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Abstract. We apply the bulk holographic dark energy in general 5D two-brane models. 
We extract the Friedmann equation on the physical brane and we show that in the general 
moving-brane case the effective 4D holographic dark energy behaves as a quintom for a 
large parameter-space area of a simple solution subclass. We find that w\ was larger 
than —1 in the past while its present value is w\ a ~ —1.05, and the phantom bound 
W\ = — 1 was crossed at z p ss 0.41, a result in agreement with observations. Such a 
behavior arises naturally, without the inclusion of special fields or potential terms, but 
a fine-tuning between the 4D Planck mass and the brane tension has to be imposed. 
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1. Introduction 

Holographic dark energy [1, 2, 3, 4] is a recently developed ingenious idea of explaining the 
observed Universe acceleration [5]. It arises from the cosmological application [6] of the 
more fundamental holographic principle [7, 8]. Although there are some objections about 
the applicability of holography to a cosmological framework [9] , holographic dark energy 
has opened new research directions, revealing the dynamical nature of vacuum energy by 
relating it to cosmological volumes. The background on which it is based, is the black 
hole thermodynamics [10, 11] and the connection between the UV cut-of of a quantum 
field theory, which is related to vacuum energy, and a suitable large distance of the theory 
[12]. This connection, which was also known from AdS/CFT correspondence, proves to 
be necessary for the applicability of quantum field theory in large distances. The reason 
is that while the entropy of a system is proportional to its volume the black hole entropy 
is proportional to its area. Therefore, the total energy of a system should not exceed the 
mass of a black hole of the same size, since in this case the system would collapse to a 
black hole violating the second law of thermodynamics. In holographic statistical physics 
terms this is equivalent to the exclusion of those degrees of freedom that would collapse. 
When this approach is applied to the Universe, the resulting vacuum energy is identified 
as holographic dark energy. 

Until now, almost all works on the subject have been formulated in the standard 4D 
framework. However, brane cosmology, according which our Universe is a brane embedded 
in a higher-dimensional spacetime [13, 14], apart from being closer to a higher-dimensional 
fundamental theory of nature, it has also great phenomenological successes [15]. In our 
recent works [16] we presented a generalized and restored holographic dark energy in the 
braneworld context. The basic argument was that in a higher-dimensional spacetime, 
it is the bulk space which is the natural framework for the cosmological application 
(concerning dark energy) of holographic principle, and not the lower-dimensional brane- 
Universe. This is obvious since it is the maximally-dimensional subspace that determines 
the properties of quantum-field or gravitational theory (such as cut-off's and vacuum 
energy), and this holds even if we consider brane cosmology as an intermediate limit of an 
even higher- dimensional fundamental theory of nature. To be more specific we recall that 
in braneworld models, where the spacetime dimension is more than 4, black holes will in 
general be D-dimensional [10, 11], no matter what their 4D effective (mirage) effects could 
be. Therefore, although holographic principle is itself applicable to arbitrary dimensions 
[7, 17] its cosmological application concerning dark energy should be considered in the 
maximal uncompactified space of the model, i.e. in the bulk. Subsequently, this bulk 
holographic dark energy gives rise to an effective 4D dark energy with "inherited" 
holographic nature, and this one is present in the (also arisen from the full dynamics) 
Friedmann equation of the brane. One can in general acquire either different or exactly 
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identical 4D behavior, comparing to that obtained in conventional 4D literature [1, 2, 3, 4]. 
However, even in the second case, the physical interpretation is radically different. 

In our previous works [16] we used, as a specific example, a general braneworld model 
with one brane. The arbitrary large extra dimension of this case imposes no restrictions 
on the application of bulk holographic dark energy, thus we recovered all the results of 4D 
literature. In the present paper we are interested in investigating the case where the bulk 
is finite, and for this purpose we use the well explored two-brane model, where the branes 
constitute the boundaries of the extra dimension [14, 18, 19, 20]. From the first moment it 
becomes clear that if the branes are steady, i.e. the extra dimension has a constant size, 
then the bulk dark energy loses its dynamical- holographic nature. The reason is that, 
applying the holographic dark energy arguments, we cannot consider an arbitrary large 
bulk black hole in this case. Although we could still use non-spherical exotic solutions 
such as black rings and black "cigars" [21], or highly rotational or/and charged black 
holes [22] , if we desire to maintain the simplicity and universality which lies in the basis of 
holographic dark energy we must remain in the aforementioned framework, i.e preventing 
the horizon of a spherical bulk black hole being larger than half the interbrane distance. It 
seems that holographic dark energy is in contradiction with brane stabilization mechanism 
[23]. 

In this work we examine braneworld models where the bulk is finite but with 
moving boundaries (branes). In this case holographic dark energy is applicable and 
the corresponding cosmological length is the interbrane distance. However, one still 
needs an additional fine-tuning assumption, arising form the time variation of the model 
parameters. We want to study the behavior of the effective 4D holographic dark energy, 
and especially its dependence on the metric scale factor. The rest of the text is organized 
as follows: In section 2 we present the holographic dark energy in the bulk and in section 
3 we apply it to a general two-brane model in 4+1 dimensions. Finally, in 4 we discuss 
the physical implications of the obtained results. 



2. Formulation of Holographic Dark Energy in a General Bulk 

In this section we display the basic results of the bulk holographic dark energy, formulated 
in [16]. The mass Mbh of a spherical and uncharged D-dimensional black hole is related 
to its Schwarzschild radius r s through [11, 24]: 

M BH = r^- 3 (V^M D ) D - 3 M D ^^-, (2.1) 

81 {—) 

where the D-dimensional Planck mass is related to the D-dimensional gravitational 
constant Go and the usual 4-dimensional Planck mass M p through: 

M D = G^\ 
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M 2 p = M£-V d _ 4 , (2.2) 

with Ve>_4 the volume of the extra-dimensional space [11]. 

If Pad is the bulk vacuum energy, then application of holographic dark energy in the 
bulk gives: 

PAD Vo\(S D ^) < r D -*{^M D ) D -*M D (2.3) 

81 {—) 

where Vol(iS D ~ 2 ) is the volume of the maximal hypersphere in a .D-dimensional spacetime, 
given from: 

Vol(<S D - 2 ) = A D r D ~\ (2.4) 



with 



D-1 
7T 2 



A D ~ 



(W 

D-2y 



D-1 



for D — 1 being even or odd respectively. Therefore, by saturating inequality (2.3) 
introducing L as a suitable large distance and c 2 as a numerical factor, the corresponding 
vacuum energy is, as usual, viewed as holographic dark energy: 

Pad = C 2 (v^M D ) B - 3 M D ^-^- IT 2 . (2.6) 

81 {—) 

As was mentioned in [16], the "suitable large distance" which was used in the definition of 
L in (2.6) should be the Hubble radius, the particle horizon, or the future event horizon 
[1, 25, 4], with the last ansatz being the most appropriate. However, in the case of a finite 
bulk with varying size it is obvious that L should be just this bulk size, with the reason 
being again the foundations of holographic dark energy which prevent the use of a length 
larger than that. On the other hand, if the future event horizon is smaller than the bulk 
size then one should use that instead of the bulk size. In this case our braneworld model 
does not "feel" the other bulk boundary and the picture is equivalent to the single-brane 
model investigated in [16]. 

Finally, let us make a comment concerning the sign of bulk holographic dark energy. 
In the original Randall-Sundrum model [11] the bulk cosmological constant should be 
negative in order to acquire the correct localization of low-energy gravity on the brane. 
Such a negativity is not a fundamental requirement and is not necessary in more complex, 
non-static models, like the one of the application of the next section. However, as was 
mentioned in [16] and generally speaking, holographic dark energy is a simple idea of 
bounding the vacuum energy from above. It would be a pity if, despite this effort, 
one could still have a negative vacuum energy unbounded from below, because then 
holographic dark energy would loose its meaning. If holography is robust then one should 
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reconsider the case of a negative bulk cosmological constant (although subspaces, such 
as branes, could still have negative tensions). Another possibility is to try to generalize 
holographic dark energy to negative values, in order to impose a negative bound. The 
subject is under investigation. 

3. Holographic Dark Energy in 5D Braneworld models with moving branes 

We are interested in applying the bulk holographic dark energy in general 5D braneworld 
models where the bulk is bounded by two branes. The inclusion of a second brane is 
probably the best way of eliminating possible "naked" metric singularities, thus allowing 
for more complex and realistic models [14, 18, 19, 20]. We consider an action of the form: 



In the first integral M5 is the 5D Planck mass, pA5 is the bulk cosmological constant which 
is identified as the bulk holographic dark energy, and R is the curvature scalar of the 5D 
bulk spacetime with metric gAB- The second term corresponds to two (3 + l)-dimensional 
branes, which constitute the boundary of the 5D space. 7 is the determinant of the 
induced 4D metric 7^ on them, Vi stand for the brane tensions and Cy. at is an arbitrary 
brane matter content [26]. 

As usual the two branes are taken parallel, y denotes the coordinate transverse to 
them and we assume isometry along 3-dimensional x slices including the branes. For the 
metric we choose the conformal gauge [18, 19]: 



This metric choice, along with the residual gauge freedom (t, y) — > (t', y') which preserves 
the 2D conformal form, allows us to "fix" the positions of the branes. Without loss of 
generality we can locate them at y — 0, 1, having in mind that their physical distance is 
encoded in the metric component B(t, y), and at a specific time it is given by [18, 19]: 



quantity that is invariant under the residual gauge freedom in our coordinates. The 
reason we prefer the metric (3.8), instead of the usual form in the literature, is that in the 
later case the brane positions are in general time-dependent and the various boundary 
conditions are significantly more complicated. Thus, our coordinates are preferable for 
numerical calculations, despite the loss of simplicity in the definitions of some quantities. 
Eventually, the physical interpretation of the results is independent of the coordinate 
choice. 



S = I d A xdy^g- (M 3 5 R - p A5 ) + £ / d 4 x^ (££f - V,) , (3.7) 



ds 2 = e 2B ^ (-dt 2 + dy 2 ) + e 2A ^dx 2 . 



(3.8) 




(3.9) 
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The non-trivial 5D Einstein equations consist of two dynamical: 



A- A" + 3 A 2 - 3 A' 2 = -|-3 e 2B p A5 (3.10) 



B-B"- 3A 2 + 3A' 2 = -—3 e 2B PA5 , (3.11) 

and two constraint equations: 

- A' A + B'A + A'B-A' = (3. 12) 

2A' 2 - A'B' + A" -A 2 -AB = - e 2B p A5 , (3.13) 

where primes and dots denote derivatives with respect to y and t respectively. It is easy 
to show that the constraints are preserved by the dynamical equations. 

We consider a brane-Universe, which as usual is identified with the brane at y — 0, 
containing a perfect fluid with equation of state p = wp (in the following we omit the 
index for the physical brane quantities and we keep the index 1 for the ones on the brane 
at y — 1). For the hidden brane we consider for simplicity just a brane tension, although 
we could also consider some matter-field content [27]. The reason we use a second brane is 
to eliminate possible "naked" metric singularities. Therefore, assuming S 1 /7a2 symmetry 
across each brane we restrict our interest only in the interbrane space. 

Integrating on a small y interval around the branes and using the boundary terms in 
the action we obtain the following junctions (Israel) conditions: 

|B ' lo = 3M eB ° (2p + 3p ~ v) (3 ' 14) 



for the physical brane, and 

[*l.= j^«*V. (3-15) 

for the hidden one, where B = B (t) and £>i = B\{t) are the values of B(y,t) at the 
branes at y — 0, 1 respectively. In the expressions above we use the following relations, 
resulting from S 1 /^ symmetry, for the jump of any function across the branes: 

[Q% = 2Q'(0 + ) , [Q'h = -2Q'(r). (3.16) 

Finally, the induced 4D metrics of the two ("fixed" -position) branes in the conformal 
gauge are simply given by 

d s 2 = -dr 2 + a 2 (r)c/x 2 , (3.17) 
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with 

dn = e Bi dt (3.18) 

ai = e M (3.19) 

the proper times and scale factors of the two branes (i = 0,1). Thus, for the Hubble 
parameter on the branes we acquire 



_ 1 da 
Hi = - — 
a ar 



e~ B 'A h (3.20) 



which is invariant under residual gauge transformations [18]. As we have mentioned, in 
the following we omit the index for the physical brane quantities. 

In order to acquire the cosmological evolution on the physical brane we proceed as 
follows: Equations (3.10)-(3.13) hold for the whole spacetime, including the branes. In 
the later case we have to use the junction conditions (3. 14), (3. 15) for the calculation of 
the first spacial derivatives. Therefore, eliminating A" form (3.10) and (3.13), and making 
use of (3.14) and of the Hubble parameter relation (3.20), we finally obtain: 

f + = ik »» + Sik I"" 2 + V, - %<, + V) + 2V\ , (3.2!) 

where H(t) is the Hubble parameter of the physical brane, with r its proper time. The p 2 
term on the right hand side of (3.21) is the usual term present in braneworld cosmology. 
Taking the the low-energy limit (p <C V) and knowing that conventional Friedmann 
equations give: 

dH 9 4ir , . 16-7T 

with pa = pA4 the 4D dark energy, it is obvious that brane evolution coincides with that 
derived from conventional 4D cosmology if we identify: 

Ml 

V = 487r-|, (3.23) 

p 

and 

1 Ml M9 

Relation (3.24) provides the (effective in this higher- dimensional model) 4D dark 
energy in terms of the bulk holographic dark energy, which according to (2.6) is given by: 

p A5 = c 2 ^M 3 5 L- 2 . (3.25) 

The holographic nature of pas is the cause of the holographic nature of pa- As we have 
already mentioned, in the two-brane model examined in the present work, the cosmological 
length L should be the interbrane distance L 5 , which is given by (3.9). Furthermore, using 
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(2.2) we eliminate the 5D Planck mass M5, in terms of the standard 4D M p , through 
Mf = Mp/L 5 , since L 5 is the volume (size) of the extra dimension. Note that the 
varying behavior of the extra-dimension size will give rise to a varying M p , i.e. varying 
4D Newton constant G4, and this is in general an inevitable consequence of moving-brane 
models [28, 3]. A detailed discussion on this point is given at the end of this section. 
Thus, we finally acquire the following form for the effective 4D holographic dark energy: 

PA =[^- 2+ 2^MlLf. (3.26) 

As we see, we have resulted in a simple holographic relation, despite the complicated 
nature of the model. 

Our goal is to reveal the dependence of pa on the physical brane scale factor a. Since 
we have recovered the standard Friedmann equation (through identifications (3.23) and 
(3.24)) and using (3.9) for the calculation of L5 we finally obtain: 

Therefore, despite the simple form of relation (3.26), complexity has reappeared in the 
non-localized nature of the above equation. Such a behavior was expected and is a result 
of the "global" properties of bulk holographic dark energy. This is a radical difference 
comparing to the single-brane case of [16] where the arbitrary large bulk allowed for a 
brane-localized equation form. In the present case, the solution of the full 5D equations is 
indispensable. Namely, we have to solve (3.10)-(3.13) under boundary conditions (3.14)- 
(3.15), imposing the conventional time evolution for p and p. Knowing A(t, y) and B(t, y) 
we calculate the interbrane distance L 5 (t) through (3.9) and then /?a(£) through (3.26). 
We calculate the physical brane scale factor a(t) using (3.19), and by eliminating t 
we obtain the questioning relation for p\(a). Finally, as usual, we identify p\(a) with 
Pa(o) ~ a~ 3( - 1+WA \ and we extract the form of w\(z), with z = ^ — 1 and ao the value of 
a at present time. 

The aforementioned procedure for the derivation of Pa (a) relation is impossible to be 
completed analytically (analytical solutions can be obtained only for stationary cases, i.e. 
with constant interbrane distance, which as we have mentioned are in contradiction with 
the holographic nature of dark energy). However, avoiding a full numerical approach, for 
the purpose of this work we examine the following solution class: 

y + Qi' 



A(t ,y) =B(t ,y) = In 



(3.28) 



<?2 

A(t ,y)=B(to,y)=q 3 [y + qi}, (3.29) 

which satisfy the constraint equations (3. 12), (3. 13) at t (and therefore at every t since 

(3. 12), (3. 13) are preserved by the equations of motion) provided that (1 + 2gi) 2 g| = 
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c 2 /(27r). Moreover, boundary conditions (3.14) and (3.15) are fulfilled imposing V = 
6Mlq 2 /q 2 1 and V x = -6M 5 3 g 2 /( ?1 + l) 2 . 

Investigating the low-energy (late-time) evolution of the aforementioned solution 
class, i.e. omitting p and p in boundary conditions which make them significantly simpler, 
we obtain interesting results. In particular, for a large area of the parameter-space we find 
a reasonable w\(z) form, with the basic requirement (necessary but not always efficient) 
being the decreasing of interbrane distance. Fortunately, a decreasing interbrane distance 
seems to have larger probability than an increasing one (since in the later case the system 
is often unstable or "naked" singularities do appear between the branes). Furthermore, 
there are many parameters present in the model and in solution class (3.28)-(3.29). These 
characteristics make it substantially easier to acquire a reasonable w\(z) form in the 
5D framework described in this work, than in conventional 4D holographic dark energy 
[1, 2, 3, 4]. In addition, we do not have to fine-tune the constant c in the definition of pA5 
in relation (3.25). In fig. 1 we depict w^(z) for q\ = —1.4, q 2 ~ 11, q% ~ 0.22 and c = 1, 
and using the unit M 5 = 1 (the specific scale does not affect the w^z) form). We observe 

-0.8 1 




-1.2 J , . , . , . , 

12 3 

Z 

Figure 1. w\ versus z for solution class (3.28)- (3.29), using q\ = —1.4, (72 ~ 11, 
53 « 0.22, c = 1 and the unit M5 = 1. The present value of wa is Wa ~ —1.05 and the 
phantom bound Wa = — 1 was crossed at z p w 0.41. 

that the effective 4D holographic dark energy behaves like a "quintom" [29] , that is % was 
larger than —1 in the past while its present value is 10 a ~ —1.05, and the phantom bound 
w\ = —1 was crossed at z v w 0.41. This result is in agreement with late acceleration of 
the Universe and dark energy constraints imposed by observations [30, 31]. Note that 
quintom behavior arises naturally in our higher-dimensional brane model, without the 
inclusion of extra fields or specific potential terms by hand. Moreover, as we mentioned 
above, we do not have to fine-tune the solution parameters, since a large area of the 
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parameter-space leads to similar behavior (the case of fig. 1 just corresponds to a good 
representative in comparison with observations). Finally, other parameter areas in the 
aforementioned solution subclass, as well as a general numerical investigation beyond this 
ansatz (a hard task due to various instabilities [18, 19]) reveal interesting but likely un- 
physical w\ behavior, such as chaotic oscillations with respect to z. However, this could 
still be the case in our Universe for larger z. The subject is under investigation. 

Let us finish this section with some comments on the dynamical nature of some 
quantities. As we have mentioned in the introduction, the main goal of the present work 
is the study of bulk holographic dark energy in a finite braneworld model. From the 
first moment we educe that if the bulk boundaries (the two branes) are steady then 
holographic dark energy becomes a constant, losing its dynamical nature. In other words, 
the concept of holographic dark energy is in contradiction with finite-bulk framework, 
unless we consider moving-brane models as a necessary outlet. Unfortunately, this choice 
leads to some undesirable consequences. Indeed, with a varying interbrane distance it 
is obvious from (2.2) that either the 5D Planck mass M 5 or its 4D counterpart M p , or 
even both, should change with time. Since M5 is a fundamental quantity of our model we 
desire to maintain it as a constant. Thus, M p and therefore the 4D Newton's gravitational 
constant G4 are the ones that reflect the dynamical nature of the bulk, and this is a 
common feature of moving-brane models [28, 3]. However, in order to acquire the correct 
effective 4D cosmological evolution we need to make the identification (3.23), a usual 
approach of braneworld cosmology. This relation "matches" the 4D Planck mass (which 
is determined by dimensional reduction of the 5D action) with the brane tension (which 
should be given by the vacuum-energy predicted by the effective 4D QFT), and brane- 
cosmologists hope to acquire its justification by a fundamental theory of nature, unknown 
up to now. Unfortunately, in the present model of moving branes this fine-tuning must 
hold at all times and this is definitely an additional assumption, independent of the rest 
formulation. Thus, the extension of holographic dark energy to finite-bulk models is still 
obscure, since one expects from a fundamental theory to justify the aforementioned eternal 
fine-tuning. 

In order to provide a clear picture of these features, in figure 2 we depict the evolution 
of the 4D Newton's constant G4 divided by its present value, in terms of z, for the same 
parameter values of fig. 1. We observe that G4 acquires its minimum value at z p ~ 0.41, 
which corresponds to the phantom-divide crossing of W\(z). This behavior of G4 is 
consistent with cosmological constraints which restrict its deviation between ±5% [32]. 
However, we mention here that not all numerical solutions satisfy these limits. Indeed, 
our numerical elaboration reveals that out of 10 3 solutions that present an acceptable 
(quintom) form for W\(z) (similar to that depicted in fig. 1), only « 20% correspond to 
an acceptable G 4 -behavior, too. 
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Figure 2. Evolution of the Newton's constant G4, divided by its present value, versus 
z for solution class (3.28)-(3.29), using q\ = —1.4, qi ~ 11, (73 ~ 0.22, c = 1 and the 
unit M5 = 1. The minimum value is obtained at z p at 0.41, which corresponds to the 
phantom- divide crossing ofw\{z). 



4. Discussion-Conclusions 

In this work we apply the bulk holographic dark energy in a general braneworld model 
with moving branes. Such a generalized bulk version of holographic dark energy is 
necessary if we desire to match the successes of brane cosmology in both theoretical and 
phenomenological-observational level, with the successful, simple, and inspired by first 
principles, notion of holographic dark energy in conventional 4D cosmology. In particular, 
as we showed in [16], the bulk space is the natural framework for the cosmological 
application, concerning dark energy, of holographic principle, since it is the maximally- 
dimensional subspace that determines the properties of quantum-field and gravitational 
theory, and the black hole formation. Subsequently, this bulk holographic dark energy 
will give rise to an effective 4D dark energy with "inherited" holographic nature, and this 
one will be present in the effective Friedmann equation of the brane. 

Applying the bulk holographic dark energy in the well investigated two-brane model, 
and using the interbrane distance (bulk size) as the cosmological length in its defini- 
tion, we deduce that the branes have to move in order for dark energy to preserve its 
holographic-dynamical nature. In this case we extract the effective Friedmann equation 
on the physical brane, which has a non-localized (on the brane) form since the interbrane 
distance must be calculated from the full 5D dynamics. This complexity is a result of 
the "global" characteristics of bulk holographic dark energy. Numerical investigation on 
a simple solution subclass reveals a quintom behavior [29] for w\{z). In particular, w\ 
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was larger than —1 in the past, it crossed the phantom divide w\ = — 1 at z p m 0.41, and 
its present value is w\ ~ —1.05. This behavior is in notable agreement with observations 
[30, 31] which give wa = — 1.02±q;]| and z v = 0.46 ± 0.13. However, although we have 
not included any special fields or specific potentials, an additional assumption has to be 
imposed in order for this holographic dark energy application to be valid. Namely, the 
fine-tuning between the brane-tension and the 4D Planck mass. Furthermore, Newton's 
constant acquires a dynamical nature too and one has to be careful in order to be consis- 
tent with cosmological constraints [32]. In conclusion, the investigation of the extension 
of holographic dark energy in finite-bulk models reveals that a reasonable u>a (-^-behavior 
is acquired relatively easily, only under a fundamental fine-tuning. Definitely, the combi- 
nation of holographic dark energy with brane cosmology is an interesting subject which 
needs further investigation. 
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